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ABSTRACT 

We report on a series of outbursts of the high mass X-ray binary XTE J1946+274 in 2010/201 1 as observed with INTEGRAL , RXTE, 
and Swift. We discuss possible mechanisms resulting in the extraordinary outburst behavior of this source. The X-ray spectra can be 
described by standard phenomenological models, enhanced by an absorption feature of unknown origin at about lOkeV and a narrow 
iron Kor fluorescence line at 6.4 keV, which are variable in flux and pulse phase. We find possible evidence for the presence of a 
cyclotron resonance scattering feature at about 25 keV at the 93% level. The presence of a strong cyclotron line at 35 keV seen in data 
from the source’s 1998 outburst and confirmed by a reanalysis of these data can be excluded. This result indicates that the cyclotron 
line feature in XTE J 1946+274 is variable between individual outbursts. 
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1. Introduction 

Due to the ~10 12 G strong magnetic field at the magnetic poles 
of many accreting neutron stars in high-mass X-ray binaries, cy- 
clotron resonance scattering features (CRSFs or cyclotron lines) 
are observable in the X-ray spectra of these sources. These lines 
originate from photons generated in the accretion column of a 
magnetized neutron star interacting with electrons in the column , 
since their motion perpendicular to the B - field is quantized into 
Landau-levels with energy differences 

“- i2k ' v (i<>y (i > 

To date CRSFs have been reported for about 20 X-ray pulsars 
(Caballero & Wilms 2011). 

The 15.8 s pulsar XTE J 1946+274 was first detected in 1998 
(Smith &Takeshima 1998; Wilson etal. 1998). It is a transient 
X-ray source with a Be-type companion (Verrecchia et al. 2002). 
This kind of main-sequence B stars shows Balmer emission lines 
caused by a quasi-Keplerian equatorial disk near the Be star (see, 
e.g., Hanuschik 1996; Slettebak 1988). Interaction between this 
disk and the orbiting neutron star can lead to violent X-ray out- 
bursts resulting in the appearance of a bright X-ray sources in 
the sky. 

The initial outburst of XTE J 1946+274 in 1998 lasted about 
three months. Heindl et al. (2001) reported on the discovery of 
a CRSF at an energy near 35keV during this outburst, which 
was followed by several fainter outbursts separated by ~80d 
(Campana et al. 1999). This separation was later established by 
Wilson et al, (2003) as half of the ~170d orbital period. The oc- 
currence of two outbursts per orbit could be related to the mis- 
alignment of the Be star’s angular momentum and the orbital 


plane of the neutron star. While we look nearly onto the pole of 
the Be star, the orbital inclination is £46° (Wilson et al. 2003). 

After a phase of nearly periodic flaring between 1998 and 
2001, XTE J1946+274 went into quiescence until 2010 June, 
when a new sequence of outbursts started (see, e.g., Muller et al. 
2010). The principal outburst behavior in 2010 was the same 
as that observed in the 1990s, with two outbursts per orbit. As 
shown in Fig. 1, however, the five outbursts seen during 2010 
are neither clearly connected to the times of periastron and apas- 
tron passages (based on the orbital ephemeris from Wilson et al. 
2003), nor to any other specific orbital phase. This new outburst 
episode ended in 201 1 June and the source was back in quies- 
cence at the time of writing. 

In this paper, we report on INTEGRAL , Swift , and RXTE 
observations of the first and the third outburst of this series in 
2010 June/July and 2010 November/December, respectively. In 
Sect. 2 we give a summary of the observations. In Sect. 3 we de- 
scribe the continuum model and study the time and pulse phase 
resolved behavior of the spectra. Summary and conclusions are 
given in Sect. 4. 

2. Observations and data reduction 

The 2010/2011 outburst series started on 2010 June 4, 
when Swift [BAT detected an increase of the X-ray flux of 
XTE J1946+274, rising up to 40mCrab (15-50 keV) within 
three days (Krimm et al. 2010). This first outburst lasted about 
60 d and reached a flux of ~140mCrab (see Fig. 1). The subse- 
quent four outbursts lasted between 30 and 50 d, each, reaching 
almost the same maximum flux level between 40 and 60 mCrab. 
The separation of the outbursts is between 60 and 90 d. A pecu- 
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Fig.l. 15-50keV Swift/BAI light curve of the 2010/2011 out- 
burst series. The times of periastron and apastron passages are 
marked by P and A, respectively. These epochs and the corre- 
sponding uncertainties were calculated using the orbital solution 
from Wilson etal, (2003). The insets provide a closer view on 
the first (Jun/Jul) and the third (Nov/Dec) outburst. The blue, red, 
and green arrows in these insets indicate the observations times 
for PCA, ISGRI, and XKT, respectively. Epochs over which data 
were summed for the time resolved spectral analysis are indi- 
cated with I to VI. 

liar behavior was observed during the third outburst where the 
primary maximum was followed by another brightening, reach- 
ing again luminosities up to ~50 mCrab (Fig. 1 , inset). 

We present data from the Proportional Counter Array 
(PCA, Jahodaetal. 2006) on board the Rossi X-ray Timing 
Explorer {RXTE, Bradtetal. 1993), the INTEGRAL Soft 
Gamma-Ray Imager (ISGRI, Lebrun etal. 2003). on board 
the INTErnational Gamma-Ray Astrophysics Laboratory 
( INTEGRAL , Winkler etal. 2003), and the X-ray Telescope 
(XKT, Burrows etal, 2005) on board the Swift Gamma-Ray 
Burst Explorer (Gehrels 2004). Data were reduced with the 
standard analysis pipelines, based on heasoft (v. 6.10 and 6.11) 
and INTEGRAL OSA v. 9.0. XTE J 1946+274 was monitored 
by RXTE regularly during the first and the third outburst. The 
source was also sporadically in the field of view of INTEGRAL 
during the first and the third outburst. Swift pointings were 
available only during the main and secondary peak of the third 
outburst. Table 1 contains a log of the observations with these 
satellites, which are also indicated in Fig. 1 . 

PCA consisted of five proportional counter units (PCUs) 
with a field of view of ~1°, sensitive between 2 and 90keV. 
Since PCU2 is known to be the best calibrated one (Jahoda et al . 
2006), only data from the top layer of this PCU are used. We ob- 
tained 2-60 keV light curves with 0.125 s resolution, spectra in 
the standard2f mode, and pulse phase resolved spectra using 
GoodXenon data. The light curves were corrected to the bary cen- 
ter of the solar system using f axbary 1 . The PCA background 
model SkyVLE was used for PCA background subtraction. Due 
to large uncertainties in the orbital parameters (Wilson etal. 
2003) no correction could be performed for the neutron star’s 
orbital motion. The lack of sufficient statistics prevented us from 
improving the existing orbital solution. Data from the High 

1 http://heasarc.gsfc.nasa.gov/lheasoft/ftools/fhelp/faxbary.txt 


Table 1. Summary of all observations used. 


ID“ 

start date 
(2010) 

MJD 

^exp 

[s] 

cts* 

[10 s ] 

e c 

PCA 

01-00 

Jun 20 

55367.12-55367.14 

1584 

2.92 

l 

01-01 

Jun 20 

55367. 1 8-55367 .20 

1408 

252 

I 

01-02 

Jun 21 

55368.69-55368.77 

4080 

7.69 

I 

01-03 

Jun 24 

55371.74-55371.78 

2880 

537 

1 

02-00 

Jun 26 

55373.50-55373.61 

5776 

10.51 

1 

02-03 

Jun 28 

55375.14-55375.17 

2800 

5.04 

I 

02-01 

Jun 28 

55375.66-55375.70 

3184 

5.62 

I 

02-02 

Jun 30 

55377.62-55377.66 

3200 

536 

I 

03-00 

Jul 03 

55380.89-55380.92 

2000 

3.06 

n 

03-01 

Jul 03 

55380.96-55380.99 

3056 

4.61 

n 

03-02 

Jul 05 

55382.92-55383.02 

6368 

8.96 

n 

03-03 

Jul 07 

55384.48-55384.59 

5872 

7.82 

ii 

04-00 

Jul 10 

55387.03-55387.13 

6144 

7.46 

n 

04-01 

Jul 11 

553 88 .73-55388 .77 

3184 

3.65 

h 

04-02 

Jul 13 

55390.95-55390.99 

3200 

337 

n 

05-00 

Jul 16 

55393.77-55393.80 

3248 

3.08 

n 

05-01 

Jul 16 

55393.83-55393.85 

1392 

1.28 

ii 

06-02 

Nov 23 

55523.19-55523.23 

2640 

235 



06-01 

Nov 24 

55524.17-55524.24 

3200 

2.65 

iii 

06-00 

Nov 25 

55525. 15-55525. 17 

1392 

1.21 

iii 

07-00 

Nov 26 

55526.54-55526.56 

2000 

1.76 

hi 

07-01 

Nov 28 

5552836-5552839 

2656 

235 

m 

07-02 

Nov 30 

55530.58-55530.62 

2928 

2.29 

IV 

08-00 

Dec 03 

55533.26-5553330 

3168 

238 

IV 

08-01 

Dec 05 

55535.28-5553531 

2432 

1.64 

V 

08-02 

Dec 07 

5553731-5553734 

2624 

1.60 

V 

ISGRI 

938 

Jun 20 

55367 .46-53368.44 

45174 

16.98 

I 

983 

Oct 31 

55500.76-55503.45 

35224 

-0.20 



988 

Nov 18 

55518.04-55518.16 

4125 

0.13 

— 

989 

Nov 18 

55518.71-55519.42 

10940 

1.83 

— 

993 

Nov 30 

55530.67-5553235 

54470 

6.94 

IV 

XRT 

001 

Nov 26 

55526 .60-55526 .66 

2475 

0.13 

IV 

002 

Nov 28 

5552839-55528.47 

2630 

0.12 

IV 

003 

Nov 30 

55530 39-5553 0 .48 

2424 

0.12 

V 

004 

Dec 02 

55532.00-55532.08 

1442 

0.07 

V 

005 

Dec 22 

55552 .49-55552 .57 

2600 

0.10 

VI 

006 

Dec 24 

55554.77-55554.98 

1838 

0.05 

VI 

007 

Dec 26 

55556.90-55556.97 

663 

0.01 

VI 

008 

Dec 28 

5555831-55558.99 

2110 

0.05 

VI 


Notes. {a) For PCA, the first column contains the number of the Obs- 
ID after 95032-12- and for XRT after 0031888. For ISGRI the revolu- 
tion number is listed . (b) Total background corrected counts. Negative 
count rates are due to uncertainties in the background determina- 
tion. (c) Epoch for data grouping, see text for details. ObsIDs without 
assignment in this column are different enough to other spectra dur- 
ing the same epoch that they cannot be combined with the other data 
or they were recorded without simultaneous low-energy measurements. 
We therefore exclude these data from our analysis. 

Energy X-ray Timing Experiment (HEXTE, Rothschild et al . 
1998) on board RXTE were excluded from our analysis since 
both HEXTE clusters were not rocking at the time of the obser- 
vations. The resulting uncertainties in the background determi- 
nation are too large for the purposes of this paper. 

The CdTe detector of INTEGRAL /] S GRI covers the energy 
range from ~18keV to 1 MeV (Lebrun etal. 2003). Thanks 
to the large field of view of INTEGRAL , XTE J 1946+274 
was detected several times during its recent outbursts between 
2010 June and 2011 April in observations pointed at Cvg X-l. 
We extracted ISGRI pulse phase averaged spectra for all obser- 
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Fig.2. Panels a and b show two example spectra of 
XTE J1946+274. Epoch I is the spectrum with the highest coun- 
trate while epoch IV is the only data set where all instruments 
are available. Histograms show the best fit continuum. Panels c-1 
show the behavior of the residuals when adding the additional 
components of the spectral model one by one. Numbers at the 
bottom right of each panel indicate the best-fit statistics, ^/dof . 
c and d: best fit using the Fermi-Dirac cutoff only, e and f: resid- 
uals after adding the lOkeV feature to the continuum model, g 
and h: residuals after adding the Galactic ridge emission to 
PCA (for epoch I the flux of this model component was fixed 
to the result obtained from epochs III and IV). i and j: resid- 
uals after adding the Fe K a line to the model. Since in the 
model fit of epoch IV the depth of the CRSF results to zero, 
panel j displays the final result for epoch IV. In panel k, the 
CRSF has been added to the data for epoch I. See also text 
for a discussion of the statistical significances of the Galactic 
ridge emission, the lOkeV feature, and the CRSF. 


vations of XTE J 1946+274 using the standard spectral extraction 
method of OSA 9 described in the INTEGRAL documentation 2 . 
We selected those observations of XTE J 1946+274 for which it 
was less than 12° off-axis. 

To cover the soft X-ray band we use data from Swift/XKT , a 
600 x 600 pixel CCD covering a field of view of 23f6 x 23f6 in 
the energy range 0.2— 10 keV. We extracted the data in Windowed 
Timing mode. For the source region we have choosen a circle 
around the source with a radius of -1'. The background was 
extracted using circles of the same radius at two off-source posi- 
tions about 2'5 away from the source. 

In the spectral analysis, for XRT we use data between 
1.5 keV and 7.0 keV. We binned the XRT and PCA data to a 
signal- to- noise ratio of ten. Due to calibration problems of PCA 
at the Xe L-edge, we exclude PCA data below 4.5 keV from 
our analysis. Discarding bins with S/N < 10 results in an up- 
per energy limit of typically 40 keV for the PCA. ISGRI covers 


2 http://www.isdc .unige .ch/intcgral/analysis#Documentation 
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Fig. 3. Results of time resolved spectroscopy. The units are 
as follows: F x in 10 -2 photons s _l cm -2 in the 7-15 keV en- 
ergy band, Afa in 10 22 cm~ 2 , Ef 0 « and E cut in keV, and A\ F e in 
10 -4 photons s" 1 cm -2 . 


the high energy band between 20keV and lOOkeV. We added a 
systematic error in quadrature to the PCA and the ISGRI spec- 
tra using canonical values of 0.5% and 2.0%, respectively (see, 
Jahoda et al. 2006, and IBIS Analysis User Manual). 


3. Spectral analysis 

We performed all fits using the Interactive Spectral 
Interpretation System (ISIS, Houck & Denicola 2000). In 
order to improve the signal to noise ratio of individual spectra 
we averaged the data over six data blocks in time taking into 
account the flux level and instrument availability. These epochs 
I— VI are defined in Table 2 and shown in Fig. 1. The spectra 
in epochs I and II cover the first outburst at the maximum and 
fading phase, respectively. Epochs III-V follow the maximum 
and the fading phase of the third outburst. The flux in the 
maximum level of this outburst (epoch III) is comparable to 
the fading phase of the first outburst (epoch II). The last set of 
observations is summarized in epoch VI. These XRT data cover 
the fading phase of the flare right after the third outburst. Before 
defining these epochs we confirmed that the spectral variability 
during these epochs is negligible. 


3.1. Spectral model 

As shown, e.g,, by Becker & Wolff (2007), the X-ray spectra of 
accretion- powered X-ray pulsars can be roughly described by a 
powerlaw with a high energy cutoff. In practical data modelling, 
this continuum has been approximated by several different con- 
tinuum models (see, e.g., Kreykenbohm et al. 2002, for a sum- 
mary). Here we describe the data using the so called Fermi-Dirac 
cutoff (FDCO, Tanaka 1986), given by 


FDCO(£) cc E~ r x 


1 + exp 



( 2 ) 


which has been successfully applied to other accreting X-ray 
pulsars such as Vela X-l (Kreykenbohm et al. 2008). This con- 
tinuum is modified by the CRSF, modeled as a line with a 
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Table 2. Results of the time resolved spectral analysis. 


e° 

start*' 

stop*’ 

Cts c 

Fi" 

F 2 a 

AV 

r 

w 


A Fc * 

WW 


I 

367 

378 

17.20 

6.91(2) 

■vgjgjgM 

2-11 

1 77+0.25 
84^ 

ips* 

2 14+3.25 

m 

I5pr 

8 .<rj! 

16.7-5 

19.441 
16 V 

10 ±4 



II 

III 

380 

524 

394 

529 

16.46 

2.93 



5.94J 
1.5 ±1.0 

46 ±20 
18 ±14 

Bill 

IV 

V 

VI 

530 

535 

552 

534 

538 

559 

1.68 

1.12 

2.54(3) 

1 . 88 ( 2 ) 

3.17(2) 

2.46(2) 

1.69(10) 

6.0 

7.25 

14 ±4 

i9-o-|;i 

16.7 

o*4!J2 

£1.5 

£ 1.0 

<30 

<27 

1.03/176 

0.78^/32 

1.21/135 


c b 

CpCA* 

ClSGRI* 

Cxrt' 

E r J 

Tg^ 

cr G J 

Wo* 

^CRSF 7 

TCRSF"* 

C'CRSF 7 

WW 

I 

II 

III 

IV 

0.943(20) 

1.00(7) 

0.97(6) 

0.97(4) 

1 

1 

1 

1 

0.97(3) 

0.88(7) 

0 . 88 ( 2 ) 

0.92(2) 

Q Q+0.4 

Q 0 + 8 $ 

[fe 

3 6 :£| 

3- 9 4 

3.7:| 

Era 

2 

2 

-380 + £* 

-170 + ^° 
1 -200 

17 Q +100 
1 ' U -100 

253 !*“! 

253 

253 

253 

0+10 

3+1 

-3 

£23 

£40 

0-65 + J;f 5 

0.65 

0.65 

0.65 

1 1 1 1 

V 

1 

1 

— 

— 


— 

— 


— 

— 

— 

— 

VI 

— 

— 

— 

1 

— 

~ 

— 


— 

— 


— 


Notes- Uncertainties and upper limits are at the 90% confidence level. Numbers without error bars were held fixed at the values listed. 

{C) Epoch for data grouping. MJD-55000. Total background corrected PCA counts between 7 and 15keV, in multiples of 10 5 . [d) Absorbed 
flux, in units of 10 ~ 2 photons s " 1 cm" 2 . F\ and F 2 cover the energy bands 10-20 keV and 7-15 keV, respectively. In epoch VI, the 10-20 keV band 
is not fully covered by the data and no flux value can be listed. (f) In units of 10 22 cm" 2 . ^ In units of keV. In units of 10^ photons s " 1 cm -2 . The 
centroid energy has been fixed to 6.4 keV, and the width to 10 _ 4 keV. In units of eV. ® Cpca, Cisqri and Cxrt are defined as the cross calibration 
and normalization constants for PCA, ISGRI, and XRT, respectively. ^ In multiples of 10" 2 . 


Gaussian optical depth profile 

exp(-r(E)), with r(E) = t C rsf X exp |-i j J • O) 

The equivalent width of this feature will be denoted with 

Wcrsf in the following. The spectra of some X-ray pulsars, in- 
cluding XTEJ 1946+274, also contain an absorption or emission 
like feature in the range 8-12 keV (Coburn 2001). The origin of 
this so-called lOkeV feature is still unclear, however, since it 
appears always at about the same energy, it is probably not re- 
lated to the magnetic field strength of the neutron star. We mod- 
elled this feature as a broad Gaussian absorption feature, as we 
also did for the CRSF (Eq. 3) with the centroid energy, Eq 9 
the width, cr G , optical depth, r G , and equivalent width, W G . 
Finally, interstellar absorption was modelled with an updated 
version of TBabs 3 , using abundances by Wilms et al. (2000) and 
cross sections by Vemer & Yakovlev (1995). 

In modelling the data, we have to take into account that the 
cross-normalization of the different instruments used is not per- 
fectly well known and the source might also slightly change in 
flux between the different observations. These effects were taken 
into account by cross calibration constants cxrt and cisgri, us- 
ing the PCA as the reference instrument. Furthermore, the PCA 
background was allowed to vary slightly in count rate. To ac- 
count for this imperfections in the modelling of the background, 
we introduced the constant Cb. 

Finally, the data modelling is affected by the proximity of 
the source to the plane of the Galaxy. Galactic ridge emis- 
sion (GRE, see, e.g., Worrall et al. 1982; Warwick et al. 1985) 
manifests itself through the presence of an emission feature at 
~6-7keV in the PCA spectrum caused by unresolved Fe K a flu- 
orescence lines. These lines are not present in the XRT data in 
epochs III and IV and therefore must be due to diffuse emis- 
sion that is picked up by the PCA only, due to its larger field 
of view. In our modelling of the PCA data we therefore intro- 
duced a model for the Galactic ridge emission based on that 
of Ebisawa etal. (2007), who described the ridge emission as 
the sum of two bremsstrahlung components and an iron line 

3 see http://pulsar.stemwarte.uiu-erlangen.de/wilms/researchAbabs/ 


complex modelled by three Gaussian lines at 6.4 keV, 6.67 keV, 
and 7.0 keV, with equivalent width ratios of 85:458:129, re- 
spectively. Since we used the PCA data down to 4.5 keV only, 
we did not account for the two soft bremsstrahlung compo- 
nents and modelled the Galactic ridge emission as the sum 
of three narrow Gaussians, with fixed energies and equiva- 
lent width ratios according to Ebisawa et al. (2007). This com- 
ponent was applied only to the PCA data and absorbed by 
the interstellar medium using the Galactic N u as determined 
from the Leiden/Argentine/Bonn (LAB) Survey of Galactic H i, 
Nh = 9.4 x 10 21 cm -2 (Kalberla et al. 2005). We determined the 
flux of the Galactic ridge emission from simultaneous fits to 
the XRT and PCA data from epochs III and IV. In epoch IV, 
residuals around 6-7 keV are only visible in PCA, while in 
XRT neither residuals from this emission, nor from a source 
intrinsic iron Ka line are detected. Thus, the residuals in the 
epoch IV PCA spectrum must be caused by Galactic ridge 
emission. In epoch III, on the other hand, we find weak ev- 
idence for the presence of a source intrinsic iron Ka line in 
the XRT spectrum. The unabsorbed flux of the 6.4 keV iron 
line of the Galactic ridge emission in both spectra is consis- 
tent with each other ((7.2 ± 2.0) X 10" 5 photons s -1 cm -2 , and 
(6.6 ± 2.0) x I0" 5 photons s -1 cm~ 2 for epoch III and IV, respec- 
tively). While these values are slightly higher than the 1.22 X 
10 -5 photons s _1 cm~ 2 reported by Ebisawa et al . (2007) , the dif- 
ference is still within the typical variation of the Galactic ridge 
emission over the Galactic plane (Yamauchi et al. 2009). Since 
the flux of the Galactic ridge emission is constant over time, 
we added this model component with parameters fixed to the 
mean value as obtained from epochs III and IV to the PCA 
spectra of all epochs. We find that for the late part of the out- 
burst, the ridge contributes 1.5% of the 3-10keV flux and 8% 
in the Fe band (6-7keV). To estimate the significance of the 
Galactic ridge emission, we performed Monte Carlo simula- 
tions of the best fit model without this feature to create a set 
of 1000 synthetic spectra. We then performed the fit allowing 
all model parameters, including the Galactic ridge emission, 
to vary. For both epochs III and IV none of these simula- 
tions led to a fake spectrum for which the improvement in^ 2 
was as large as in the real data, i.e., the probability that the 
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10 20 30 

energy [keV] 


Fig- 4. a: PCA spectra of time bins I-V (shifted in y-direction 
for better visibility), b: Ratio of the background corrected time 
resolved spectra. 

Galactic ridge emission is real is greater than 99 S % (>3.3o-). 
The remaining residuals at 6.4 keV can be explained by a nar- 
row source intrinsic iron Kor flourescence line. We model this 
feature by a thin Gaussian emission line with fixed centroid 
energy = 6.4 keV and width <XF e = 10 -4 keV. The flux, 
Ap e , (and thus equivalent width, Wpe) were allowed to vary. 
In summary, die model used can be written as 

M = TBabs X (FDCO + Fe^^ j- e v) x Gio kev x Gcrsf + GRE. (4) 


3.2. Time resolved spectroscopy 

In this Section we describe the time resolved behavior of the 
spectral parameters. For each of the six epochs we fitted the re- 
spective spectra separately, including the Galactic ridge emis- 
sion as a constant component as discussed above. Example spec- 
tra of two epochs together with the best model fit are shown in 
Fig. 2. The free fit parameters are summarized in Table 2 and 
displayed in Fig. 3. We calculated fluxes in the energy band 
10-20 keV for the epochs I-V. In order to be able to compare 
the source fluxes of all epochs, but avoid excessive extrapola- 
tion of the models, we also derived the fluxes in the energy band 
7-15 keV. We excluded the contribution from the Galactic ridge 
emission from this flux. 

The resulting for all fits does not exceed 1.3 (see 
Table 2). For epoch I, the value of x$ ed is rather high. Here, 
the greatest contribution to originates from residuals of the 
ISGRI data which are caused by calibration problems, so we ac- 
cept this fit. Note that not all model components are necessary to 
describe the data in the spectra with low statistics, e.g., epoch V 
provides an statistically too Iow;^ d of 0.78. The reason for in- 
cluding these components in these fits as well, even though the 
components overdetermine the fit model, is that this way it can 
be shown that these spectra are at least consistent with the full 



pulse phase ip 


Fig. 5. PCA-Pulse profile of epoch I (in the full PCA energy 
band, i.e., 2-60 keV), shown twice for clarity. The count rate in 
each bin is normalized to its variance relative to the mean pulse 
profile count rate. The phase bins are marked as follows. MP: 
main peak, SP: secondary peak, off: off state. 


model. In addition, we note that an overestimation of the sys- 
tematic error would also yield a^ re d that is too low. In epoch VI 
only XRT data are available and the continuum parameters are 
badly constrained from these data alone. We therefore fixed Ef 0 id 
and Evxi to the value obtained from epoch I because a change of 
these parameters affects mainly energies not covered by XRT. 
However, possible influences of these fixed parameters to the 
free fit parameters Ah and T might affect the results. 

The behavior of the photon index, T, the cutoff energy, 
Ecut, and the folding energy, £f 0 id, yields information about 
the evolution of the spectral continuum. In most cases there 
are no or only slight variations of these parameters apparent. 
However, we know from previous observations that these param- 
eters can also be significantly correlated to each other. We there- 
fore also derived a model independent illustration of the spec- 
tral changes by dividing the background subtracted PCA spectra 
from epochs II-V by the spectrum from epoch I (see Fig. 4) . The 
ratios II/I, III/I, and IV/I appear to be mainly constants, meaning 
that the variations of the continuum parameters in epochs I-IV 
are probably artificial and due to cross correlations. In contrast 
to this, ratio V/I shows a spectral softening for epoch V, caused 
by a real change of the continuum parameters. 

Except for epochs covered by XRT, the hydrogen column 
density Ah is not well determined and the best-fit parameters are 
consistent with constant Ah- In addition to the line caused by 
Galactic ridge emission, an Fe Kor fluorescence line at 6.4 keV 
is required during the first outburst and the maximum phase of 
the third outburst (epochs I— III) . The flux of this line is correlated 
with the X-ray flux of the source, as expected for a fluorescent 
line. Furthermore, for these epochs at highest luminosities 
(1-ni), the equivalent width W&, stays roughly constant. In 
the fading phase of the third outburst and its subsequent flare 
(epochs IV-VI) , this feature is consistent with zero, meaning that 
the additional emission at these energies can be explained by 
the Galactic ridge emission. To estimate the significance of the 
Fe Ka line, we performed similar Monte Carlo simulations 
as those done for the Galactic ridge emission. For epochs I 
and II, the probability that there is a source intrinsic Fe K a 
line is greater than 99.9 % (>3. 3 &). During epochs III and IV, 
where the source was much fainter, these simulations show 
that the probability for Fe line emission from the source is 
985% (2.4a*) and 62 . 0 % (0.9<x), i.e., here the Fe line region 
is dominated by Galactic ridge emission. 
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Fig. 6. a: PCA spectra of the three phase bins together with the 
pulse phase averaged spectrum (shifted in y-direction for better 
visibility), b: Ratio of the background corrected pulse phase re- 
solved spectra. 


We find no evidence for the centroid energy, the width, 
and the optical depth of the 10 keV feature to be variable 
over time. Due to the lack of statistics, the width of the fea- 
ture, <ro, in epoch III and IV cannot be constrained and is there- 
fore fixed to the mean value obtained from epochs I and II, i.e., 
<tg = 2.0 keV. In epoch V, also due to the lack of statistics, this 
feature is not required to describe the data. We therefore omit it 
from our model in this epoch, and also from epoch VI, where 
our coverage exists only below 7 keV. 

One of the most interesting questions is whether there is a cy- 
clotron line present as in the 1998 outburst (Heindl et al. 2001). 
We find a possible CRSF at -25 keV during epochs I-IV, i.e., 
those epochs where good coverage exists above lOkeV. The 
statistics during the maximum of the first outburst (epoch I) 
are good enough to obtain CRSF parameters that are well con- 
strained. As is fairly common for CRSF fits (Coburn 2001), the 
width and depth of the line are strongly correlated. We therefore 
set a lower limit of 05 keV for the width, which is comparable 
to the resolution of the PCA at these energies. Monte Carlo sim- 
ulations as described above lead to a 93% (1.81cr) probability 
that the CRSF found in epoch I is real. 

The dependence of the CRSF’s parameters to the choice of 
the approach of modeling the continuum is an important point 
which also has to be discussed. The results presented here were 
obtained by first adding the lOkeV feature to the model, and 
then accounting for the CRSF. Doing this vice versa for epoch I, 
we find another minimum of , leading to different cyclotron 
line parameters, i.e., £crsf = 29.6+ j;* keV, t C rsf = 0.34^;^, 
and ctcrsf = 6.8+^keV. The main problem for this alternative 
approach is the need of fixing parameters to certain values for 
the final fit. For example, allowing ct C rsf for the final fit to vary 
leads to an unrealistically broad CRSF which effectively models 
part of the exponential rollover and not the line. Furthermore, 
the quality of this fit is slightly worse (x 2, « 97 vs.;^ 2 « 86). We 
note that there is a third solution in which the centroid energy 
of the CRSF is in agreement with Heindl et al. (2001). However, 


Table 3. Results of the pulse phase resolved spectral analysis. 


parameter MP 



0.50-0.91 

r 

0.87 ± 0.05 

£foid [keV] 

6.7 ±0.4 

E c ut [keV] 

19.2 + j 2 

a Fc no- r 

11.9 ±2.5 

Wp e [eV] 
E g [keV] 

** 
' *-0.7 

To 

0-034!°;“* 

Wg [eV] 

-200^° 

Tcrsf 

0.12 ±0.08 

Wc* SF [eV] 

-180 ±90 

*y d ° f 

0.92/56 


SP 


0.13-034 

049!™ 

n q+04 
' ' - 0.5 
11* 

8*7 3! 

46 ±15 
9.5 ±0.3 
0 . 080 ^ 
-not! 
o.04$g 


1.10/55 


off 

0.34—0,50 & 0.91-0.13 
0 71+0.07 

7.8 ± 0.5 
15.0 + f^ 

93 ±2.2 
54 ± 14 

9.9 ± 0.2 

0.1 14 ±0.018 
-600^ 

0.13 ±0.09 
-210 ±130 
0.93/56 


Notes, Uncertainties and upper limits correspond to the 90% confi- 
dence level for one parameter of interest. 

(a) In units of photons s -1 cm -2 . The centroid energy has been fixed to 
6.4 keV, and the width to lO^keV. 


contrary to Heindl et al. (2001), this third solution has an un- 
physically broad and shallow shape and is thus not physically 
meningful. All results presented in this paper are based on the 
approach first adding the lOkeV feature, and then the CRSF, 
Finally, we note that there are two further solutions, where 
the CRSF is located at 30keV and 40keV and a line width 
of less than 1 keV. These two solutions are indeed physically 
meaningful, however, a Monte Carlo estimation of the signif- 
icances of these features yields probabilities of 85% (1.44<r) 
and 86% (1.48<r), respectively, much lower than for the so- 
lution with the CRSF at 25 keV. Furthermore, these two so- 
lutions are. supported only by about three data bins per in- 
strument (the 40 keV solution even only by ISGR1), while the 
20 keV solution is based on a much larger number of spectral 
bins. 

Including a CRSF in the fainter observations in epochs II- 
IV does not significantly improve the fits. All observations are 
in principle consistent with the presence of a weak line as that 
seen in epoch I. Including for consistency such a feature in the 
spectral modeling, fixing the centroid energy and width and leav- 
ing the optical depth respectively the equivalent width of the 
line as a free parameter effectively gives an upper limit for the 
depth of the line in these observations. Not unexpectedly, the 
limit for tcrsf becomes larger for the fainter phases of the 
outburst (Thble 2). We find no evidence for the equivalent 
width Wcrsf to be variable over time. 

3.3. Pulse phase resolved spectroscopy 

Spectral parameters are known to be variable as a function 
of pulse phase, e.g., the parameters of cyclotron lines (e.g., 
1A 1118-61; Suchy etal. 2011). To investigate such changes, 
we performed pulse phase resolved spectroscopy on epoch I 
with PCA, based on a pulse period ephemeris found by epoch 
folding. Epoch I provides the best statistics for this analy- 
sis. The ephemeris can be described with a constant spin 
up, P pulse = (-3.0 ± 0.3) x 10“ 9 ss -1 with pulse period between 
15.755s and 15.767s. The uncertainties of the pulse periods 
were typically less than 0.003 s. Since it was not possible to cor- 
rect the light curve for the orbital motion of the system, this spin 
up trend might be caused by orbital effects and not by transfer 
of angular momentum of the accreted matter onto the neutron 
star. We extracted pulse phase resolved spectra using individual 
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Fig. 7. Flux of the iron K a fluorescence line against the 
7—15 keV flux. The crosses correspond to the data from the 20 10 
outburst, and the diamond to the average spectrum of the 1998 
outburst , obtained by reextracting and analyzing the old PC A 
and HEXTE data from 1998. 


pulse periods following this ephemeris. The resulting pulse pro- 
file (Fig. 5) shows a double peaked structure, typical for such 
sources (see, e.g., Bildsten et al. 1997). 

As in the time resolved case, we used the FDCO model to fit 
the phase resolved spectra using three phase bins as indicated in 
Fig. 5: The main peak (MP), the secondary peak (SP) , and the off 
state (off). The pulse phases of the bins used are summarized in 
Table 3. Constant Galactic ridge emission is taken into account 
in the same way as for the pulse averaged analysis. The back- 
ground scaling factor for PCA Cb is assumed to be equal to the 
result from the pulse averaged analysis and was therefore frozen 
to this value. Initial fits of the pulse phase resolved spectra lead 
to large unconstraints in Nh and the width of the lOkeV feature 
cr G . These values were therefore frozen to the results form the 
phase averaged analysis. The signal to noise ratio in the spectra 
from the individual pulse phases is too low to allow us to study 
the behavior of the CRSF. Similar to the flux dependent analysis, 
*ve freeze the CRSF centroid energy and width to the result from 
the pulse averaged analysis and allow only the optical depth re- 
spectively the equivalent width to vary (due to the marginal 
significance of the CRSF no further attempt was made to study 
a phase dependence of the cyclotron line). 

Figure 6 shows the ratios of background corrected pulse 
phase resolved spectra. The ratio MP/ofF appears mainly as a 
constant, whiie the other two ratios reveal more complex shapes. 
This constant ratio is in agreement with the spectral parameters 
of the main peak and the off state . These parameters show almost 
no or only slight differences to each other. The ratio MP/SP, a 
falling line, can be explained by a change in photon index. This 
explanation is confirmed from the spectral analysis. While £ cut 
and £f ol d show only small differences, the photon index T is sig- 
nificantly lower in the secondary peak, i.e., the spectrum in this 
phase bin is harder than in the other phase bins. The centroid en- 
ergy of the lOkeV feature does not vary with pulse phase while 



i . i i i i i i i_d 

20 40 60 80 100 

-Ecrsf [keV] 


Fig. 8. a: Difference of x 2, of the best fit of epoch I with frozen 
CRSF parameters (Fig. 2k) and the x 2, of the best fit without 
CRSF. The centroid energy has been varied in steps of 1 keV. 
The red curve displays the x 1 of the best fit when using tcrsf 
and o~ crsf from Heindl etal. (2001). For the light blue curve 
the parameters from this work have been used. The dark blue 
curve shows this difference when both, the width and depth 
(the upper limit has been set to 3) of the CRSF are allowed to 
vary, i.e., the same situation as that of Fig. 2. The two addi- 
tional minima inx 1 are discussed in Sect. 3.2. The horizontal 
dashed line indicates Ax 1 = 0, and the horizontal dotted line 
seperates the logarithmic and the linear scaling in y-direction. 
The vertical dashed lines shows the position of the CRSF de- 
rived in Heindl et al. (2001) and this work, b and c: results for 
the optical depth and the width of the CRSF, respectively, 
when allowing these parameters to vary. The plateau of the 
curve displayed in panel b at ~60-70keV is caused by the 
upper limit for tcrsf of 3. 


its depth weakens during the main peak. The equivalent width 
remains constant in the main and the secondary peak, while 
it shows a larger value in the off state. Finally, we note that 
there is neither evidence for changes of the Fe line flux nor for 
the equivalent width over the pulse phase. 


4. Results and conclusions 

4.1. Outburst series 

In this paper we analyzed quasi-simultaneous XRT, ISGRI, and 
PCA observations of two outbursts of XTE J1 946+274 during a 
series in the second half of 2010. 
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Before the onset of this outburst series, XTE J 1 946+274 was 
in a state of quiescence for almost one decade, i.e., about 20 or- 
bits of the neutron star. Even though the formation, structure, 
and dynamics of Be disks are even today far from being com- 
pletely understood (see, e.g.. Draper et al. 201 1 , and references 
therein), the missing mass accretion onto the neutron star during 
that time is probably due to the absence of a sufficiently large Be 
disk during that time. A new outburst was only possible once the 
disk had been replenished and accretion could be triggered. 

This series shows similar behavior as the one observed in 
1998, e.g., two outbursts are observed per orbital period. The 
outbursts of XTE J1 946+274 could also be similar to those seen 
in GX 301-2, which have been extensively modeled by Leahy 
(2002). These authors posit that an additionally stream of matter 
is flowing from the primary, and that a second outburst per orbit 
could be caused by the passage of the neutron star through this 
stream. These flux peaks of GX 301-2 occur near the apsides of 
this system. 

However, the outbursts of XTE J 1 946+274 do not clearly co- 
incide with the times of periastron and apastron passages of the 
neutron star (Fig. 1). This could be explained by 3-dimensional 
simulations, which show that the disturbance of the Be disk by 
the gravitational field of the neutron star could lead to a strong 
asymmetric structure of the circumstellar material, which could 
also lead to multiple X-ray outbursts during one orbital period 
( see , e .g . , Okazaki et al . 20 1 1) . 

Alternatively, the outbursts can be triggered by the neutron 
star passing through the Be disk due to a misalignment of the 
orbit and the Be star’s equatorial plane (Wilson et al. 2003). 

Another possible explanation for the irregularity of the out- 
bursts is that they are triggered by density variations in the Be 
disk and not by orbital effects alone. This assumption could be 
verified by regular optical monitoring of the Balmer Ha line, 
which is an indicator for the presence of such a disk (see, e.g., 
Grundstrom et al. 2007, and references therein). 

42. Spectroscopic results 

An absorbed Fermi-Dirac cutoff powerlaw together with an iron 
K a fluorescence line, an iron line complex between 6 and 7 keV 
caused by Galactic ridge emission, a Gaussian like absorption 
feature around lOkeV, and a cyclotron line at ~25keV repro- 
duce the observations well in terms of . We find time as 
well as pulse phase dependent variability of the continuum pa- 
rameters of XTE J1946+274. In the time resolved case, these 
changes might be caused by different accreting mechanisms de- 
pending on the mass transfer onto the neutron star. The periodi- 
cally changing line of sight with respect to the X-ray emitting re- 
gion, caused by the rotation of the neutron star, is likely respon- 
sible for the variabilities observed in the pulse phase resolved 
analysis. Furthermore, this rotation together with differences in 
the accretion geometries at the two magnetic poles of the neu- 
tron star might lead to the asymmetric pulse profile of the X-ray 
pulsar (see Fig. 5). 

Studying the behavior of the spectral shape is problematic, 
because the spectral parameters photon index T, folding energy 
£r old , as well as the cutoff energy £ cul , which describe the broad- 
band shape of the X-ray spectrum, show strong cross correla- 
tions to each other. To account for this, we calculated ratios of 
the spectra for both, the time resolved and the pulse phase re- 
solved studies. In the time resolved case, we find clear devi- 
ations from a constant ratio for epoch V. During epochs I— IV, 
the overall spectral shape seems to remain relatively constant. In 
the pulse phase resolved analysis, the spectrum of the secondary 


peak turns out to be significant harder than these from the main 
peak and the off state. Such properties of X-ray spectra are due 
to, e.g., the temperature of the visible part of the X-ray emitting 
region, in particular the accretion column. 

Another feature required to get a good fit is a source intrin- 
sic Gaussian iron Ka emission line at 6.4 keV. The disagree- 
ment in the strength of the required emission in this energy 
range between PCA and XRT in initial fits can be solved by in- 
ducing the emission of the galactic ridge. The strength of the 
Galactic ridge emission is about four times larger than found 
in other regions of the Galaxy (Ebisawaet al. 2007). Such a 
difference is consistent, however, with the typical spatial varia- 
tions of the ridge emission (Yamauchi et al. 2009). Furthermore, 
Kiihnel et al. (2012) found an emission strength consistent with 
our result for GRO J 1008-57. 

After taking into account the contribution due to Galactic 
ridge emission, the flux of the source intrinsic Fe Ka emis- 
sion line is correlated with the X-ray flux Fx (see Fig. 7), 
as also observed in other X-ray transients (see, e.g., Inoue 
1985). The Fe Ka line is significantly detected at source fluxes 
(7-15 keV) greater than about 0.035 photons s" 1 cm" 2 . At fluxes 
below this value, the source intrinsic line is consistent with zero. 
The equivalent width stays constant for high luminosi- 
ties, i.e., epochs I-Ul. Reanalyzing the earlier RXTE data 
(Heindl et al. 2001) with the continuum model employed 
here yields an equivalent width which is also in agreement 
with the result from the 1998 outburst (W^ 998 = 59 ± lOeV). 
The different values for W? c for the fainter observations, 
where the Fe Ka line is only marginally detected, can be 
caused by the uncertainties of the Galactic ridge emission, 
which significantly contributes to the data at these energies. 
In the pulse phase resolved analysis we find no evidence for vari- 
ability of the line flux. This result indicates that it is emitted in 
a region large compared to the distance travelled by light during 
one pulse period. The equivalent width also shows no varia- 
tion, which we would actually expect for a constant model 
component at varying flux levels. However, this can be ex- 
plained by the quite large relative uncertainty of Wr, which 
is on the same order of magnitude as the respective varia- 
tions. 

The lOkeV feature is present during all epochs and pulse 
phases. It shows relatively constant results for line energy, width, 
and depth in the time resolved case. The reanalysis of the 
1998 data yields an equivalent width of W^ 998 = -250^ 
which also equals the results for the current outburst. While 
we find almost no connection of the parameters of this fea- 
ture with the luminosity, its optical depth and equivalent 
width varies with pulse phase. Wq is constant daring the 
main and the secondary peak, but it strongly increases dur- 
ing the off state. This behavior indicates that the fractional 
amount of absorbed flux related to the continuum level is 
constant for the main and the secondary peak and changes 
in the off state. This behavior of the lOkeV feature could 
give rise to speculations about possible physical processes 
as, e.g., pulse phase dependent absorbing processes, produc- 
ing such a feature. On the other hand, the residuals around 
lOkeV could also be due to the failing of the spectral broad- 
band continuum models resulting in a wrong description of 
the data around these energies. More quantitative analyses 
of this feature, also from other sources where it occours in 
emission, are urgently needed to reveal the true nature of 
this enigmatic feature. 
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4.3. Cyclotron resonance scattering feature 

We find weak evidence for the presence of the CRSF, first dis- 
co\ered by Heindl et al. (2001) during the 1998 outburst series. 
The cyclotron line improved the model fit for the high signal to 
r.oise epoch I spectrum, where PCA and ISGRI data are avail- 
able. Later spectra are consistent with the presence of a CRSF 
with unchanged parameters, however, the line is not formally 
detected in these observations because of their lower signal to 
noise ratio. We stress that due to its relative weakness, the inclu- 
sion of the CRSF in our fits does not affect our results for the 
continuum parameters or for the pulse phase resolved analysis. 

If the identification of the line is correct, then its depth 
h 2010 was significantly lower than that measured in the 
1998 outburst (t*”° = 0.09!™ vs. r™ = 0.33^, 

Heindl etal. 2001). Furthermore the centroid energy of the 
CRSF in 2010 (E™, = 25.3^ keV) is significantly lower 
than in 1998 = 36.2^ keV, Heindl etal. 2001). 

Reanalyzing the earlier RXTE data as described above 
yields CRSF parameters which are consistent with those 
found by Heindl et al. ( 2001 ). Furthermore, the equivalent 
width in 1998 was also significantly larger than in the cur- 
rent data (W**® = (-2.1™) x 10 3 eV). 

One reason why constraining the CRSF is so much eas- 
ier in the 1998 data is due to the fact that, while the expo- 
sure times of the 1998 and the 2010 data are comparable, the 
10-20 keV X-ray flux of the 1998 observation is -40% higher 
than in the 2010 observation (-9.6 x 10 -2 photons s _1 cm -2 vs. 
-6.9 x 10 -2 photons s -1 cm -2 ). However, fits to the 2010 data in 
which we fix the CRSF parameters to their 1998 values do not re- 
sult in a satisfactory description of the data = 1.55). Could 
this indicate that the energy of the line varied? Assuming that 
the CRSF width and depth remained the same, we searched for 
a line similar to that found in the 1998 data. Varying the CRSF 
energy in steps of a 1 keV and refitting the continuum parame- 
ters did not yield any satisfactory results, with best-fit >1.5 
over the 20 to 60keV band (see Fig. 8, red curve). This result 
iidicates that irrespective of the difference in signal to noise be- 
tween both data sets the CRSF must have varied between both 
outbursts, and a strong line as that seen in 1998 is not consistent 
with the data analyzed here. 

Taking the 2010 CRSF values at face value, a possible ex- 
planation for the difference between the 1998 and 2010 out- 
bursts could be the flux dependence that is seen in some other 
CRSF sources (see, e.g., Caballero & Wilms 2011, for a re- 
cent discussion) . XTE J 1946+274 could therefore be a cyclotron 
source with an overall positive correlation between the X-ray 
flux and the CRSF energy, similar to Her X-l (Staubert et al . 
2007). In models for the change in cyclotron line energy in 
neutron stars this change is generally interpreted as a change 
in height of the line producing region as the accretion column 
changes with mass accretion rate (and thus luminosity). Using 
the dipole approximation, and the assumption that in 1998 the 
CRSF was emitted from the neutron star's surface (rns = 10 km 
and /tins = 1.4 Mo) yields a height difference of 1 .7 km for the 
regions where the CRSF is generated. This result is in agree- 
ment with typical estimates of several kilometers for the height 
of the accretion columns (see, e.g., Basko & Sunyaev 1976). 
Becker etal. (2012) show that a positive correlation between 
flux and energy is possible in the luminosity range where the 
slopping in the accretion column is dominated by Coulomb 
braking. For cyclotron line energies around 30keV Coulomb 
braking is the dominating braking process in the luminosity 
nnge 1-5 X 10 37 erg s -1 . For higher luminosities, radiation brak- 


ing dominates. XTE J1 946+274 *s peak outburst luminosity of 
4.5 X lO^ergs -1 is barely consistent with this range, however, 
given that the source seems to be located at the transition be- 
tween Coulomb and radiation braking and was brighter in 1998, 
the very large change in CRSF energy seems unlikely to be due 
to a pure mass accretion rate effect. 

In conclusion, while the possibility of a luminosity depen- 
dent CRSF is intriguing, the poor signal to noise ratio of the 
2010 data does not allow a definitive answer concerning the lu- 
minosity dependence of the CRSF. Further, longer, monitoring 
observ ations of XTE J 1946+274 during its next outburst episode 
are urgently needed to resolve this question. 
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